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THE APPLICATION OF GREEN’S THEOREM TO THE SOLUTION 
OF BOUNDARY-VALUE PROBLEMS IN LINEARIZED 
SUPERSONIC WING THEORY 1 

By Max. A. Heaslet and Harvard Lomax 


SUMMARY 

Following the introduction of the linearized ‘partial differen- 
tial equation for nonsteady three-dimensional compressible 
flow, general methods of solution are given for the two- and 
three-dimensional steady-state and two-dimensional unsteady- 
state equations. It is also pointed out that, in the absence of 
thickness effects, linear theory yields solutions consistent with 
the assumptions made when applied to lifting-surface problems 
for swept-back plan forms at sonic speeds. The solutions of the 
particular equations are determined in all cases by means of 
Green’s theorem and thus depend on the use of Green’s equivalent 
layer of sources, sinks, and doublets. Improper integrals in 
the supersonic theory are treated by means of Hadamard’s 
“ finite part” technique. 

Four applications of the general solutions are given: First, 
the angle-of-attack load distribution for a supersonic, yawed, 
triangular plate with subsonic leading edges is determined. 
Second, downwash is calculated along the center line in the 
plane of the unyawed triangular wing. Third, the growth of 
load distribution is presented for subsonic and supersonic two- 
dimensional flat plates either starting from rest at a uniform 
velocity or experiencing an abrupt angle-of-attack change. The 
transient effects on lift-curve slope are then calculated. Finally 
the load distribution and lift-curve slope of a specific swept- 
back lifting surface are determined at afree-stream Mach number 
of one. 

INTRODUCTION 

If the effects of viscosity are assumed small and shock-free 
compressible flow is considered, the velocity field about a 
two- or three-dimensional body placed in a uniform free 
stream is irrotational and thus possesses a velocity potential. 
In the determination of the pressures exerted on such a body 
or in the calculation of the induced velocity components, the 
theoretical aerodynamicist is concerned essentially with 
finding the velocity potential of the flow field and, thus, 
must determine the solution of a second-order nonlinear 
partial differential equation subject to certain boundary 
conditions. The known mathematical difficulties that arise 
in the treatment of such a problem make it expedient to 
resort to simplifying assumptions. In applied aerodynamics, 
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however, efficiency of flight at high speeds has focused atten- 
tion on bodies inducing relatively small velocities throughout 
the field of flow and, as a consequence, the demands of engi- 
neering furnish a guide for the mathematical simplification 
of the theory. The so-called linearized theory of compressi- 
ble flow was developed to solve such problems and, although 
considerable work of a more precise nature has been presented 
in two dimensions, a large amount of investigation in 
unsteady flight and in three-dimensional wing theory remains 
to be completed within the framework of the simplifying 
conditions. 

The present paper is restricted to a discussion of wing 
theory subject to the assumptions of linearized compressible 
flow. It therefore employs solutions of Laplace’s equation 
and the wave equation for cases where the boundary condi- 
tions are specified in the plane of the wing. Attention will 
be directed primarily to the analysis of steady-state condi- 
tions although an equivalence will be established between the 
two- dim ensional differential equation containing the time 
variable and the equation applying to three-dimensional 
supersonic wing theory. Solutions in all cases will be 
obtained through the use of Green's theorem and the resultant 
concept of Green's equivalent layer of sources, sinks, and 
doublets. The correspondence between the theoretical devel- 
opment for subsonic and supersonic speeds is particularly 
useful since experience related to analysis in either flight 
regime is more readily transferred- 

In view of the widespread use of sources, sinks, and doub- 
lets in low-speed studies and the fact that the earlier applica- 
tions to supersonic wing theory by Prandtl (reference 1) 
and Schliehting (reference 2) corresponded to the use of 
Green's equivalent layer, it is notable that later emphasis 
shifted to other methods of solution. Sources alone were 
used by Puckett (reference 3) to create symmetrical non- 
lifting wings and were also applied to the study of lifting 
triangular wings with supersonic leading edges, hut the use 
of source, sink, and doublet sheets has not been as extensive 
as might have been expected. This anomaly is even more 
apparent in view of the vast mathematical and physical 
literature centering around the use of Green’s theorem. One 
possible explanation may stem from the fact that the interest 
of the mathematician and physicist in the wave equation 
has arisen in connection with problems in acoustics, optics, 
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and vibrating membranes. Such problems introduce bound- 
ary conditions of the Cauchy type, that is, initial conditions 
need to be known both for the unknown function and its 
rate of change. The supporting surface for such boundary 
conditions cuts the characteristic cone of an arbitrary point 
in a closed curve and has been called by Hadamard (reference 
4) a duly inclined surface. In aerodynamics the supporting 
surface is nonduly inclined and cuts the characteristic surface 
or Mach cone along the arc of a hyperbola and, as a result, 
the problem is no longer of the Cauchy type and the analysis 
becomes similar to that used in subsonic theory in the solu- 
tion of Laplace's equation. Prior to the interest of the theo- 
retical aerodynamicist in supersonic wing theory, it appears 
that little attention in application was paid to this type of 
solution. 

The material presented here is divided into two main, 
divisions: Analysis and Applications. In the first part of 
the Analysis division, the linearized differential equation for 
nonsteady compressible flow is given together with the 
underlying assumptions made. Specific forms of this equa- 
tion for two- and three-dimensional steady states and two- 
dimensional unsteady states are then considered. It is 
pointed out in particular that for swept-back lifting surfaces 
linearized theory yields consistent solutions at a free-stream 
Mach number of one although the analysis of arbitrary 
thickness distributions is not possible. Following the var- 
ious equations, Green’s theorem is applied to find, in terms 
of the known boundary conditions, the desired solution by 
means of source and doublet distributions. 

Applications of the general methods are confined to four 
problems. As an example of the manner in which angle-of- 
attack load distributions are determined for a lifting flat 
plate, the case of a yawed triangular wing with subsonic 
leading edges is solved. Doublet distributions are then 
applied in the second problem to the calculation of down- 
wash behind the same wing in an unyawed position. Third, 
the growth of load distribution with time is derived for a 
supersonic two-dimensional flat plate either experiencing a 
sudden sinking motion or starting from rest at a uniform 
velocity. Such distributions are of value in the calculation 
of indicial lift functions and can be. used, together with 
Duhamel’s integral, in the study of certain dynamic maneu- 
vers. The final application considers at a Mach number 
equal to one the case of a swept-back lifting surface with 
tips normal to the free-stream direction. , 

LIST OF IMPORTANT SYMBOLS 
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speed of sound in free stream 
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lateral distance to inboard tip of swept wing 
(See fig. 9) 
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pressure coefficient 



complete elliptic integral of the second kind with 
modulus k 

complete elliptic integrals of the second kind 
with moduli — G 2 , -yjl—8 0 2 , respectively 
complete elliptic integrals of the second kind 
with moduli k h k 2i respectively 
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X—Co 


X — Cp 
Cq8q 

complete elliptic integral of the first kind with 
modulus k , 

complete elliptic integrals of the first kind with 
moduli ki, k 2 , respectively 
tan S 0 

free-stream Mach number 
local static pressure 
free-stream static pressure 


free-stream dynamic pressure 



V(3-£) 2 +/3 2 [(3/->?) 2 +(2-r) 2 ] 
V(*-0*-0*[(y-u>*- k(2-f) 2 ] 

distance traveled in half-chords 

area of wing 

time 

perturbation velocity , components parallel to 

x, y, z axes, respectively 

jump in value of u, w at the 2=0 plane 

velocity of free stream 

Cartesian coordinates 

coordinates introduced in equation (2) 

angle of attack in radians 

V|i-M 0 2 ! 


load coefficient (Gp—Cj, u ) 

semivertex angle of yawed triangular w'ing 
angle between lifting element and x axis 
angles between leading edges of yawed triangle 
and x axis 

(3 tan 5, /3 tan 5 0 , |3 tan Si 
sideslip angle of yawed triangle 

M ach angle (^arc sin 

Cartesian coordinates 
free-stream density 
region of integration in equation (10) 
perturbation velocity potential 



sign denoting “finite part” of integral 
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SUBSCRIPTS 

u subscript denoting value of variable on upper surface of 
wing 

l subscript denoting value of variable on lower surface of 
wing 

ANALYSIS 

THE PARTIAL DIFFERENTIAL EQUATIONS 

Basic differential equation. — Consider an aerodynamic 
body flying at an arbitrary Mack number M a in air initially 
at rest. If a Cartesian coordinate system x,y,z is fixed rela- 
tive to tke body, tke body may then be assumed stationary 
and situated in a free stream with tke same Mack number. 
If tke free-stream velocity vector is parallel to and in tke 
direction of tke positive x axis and if 4> denotes tke perturba- 
tion velocity potential for isentropic flow, tke linearized 
partial differential equation for $ may be written in tke form 

(Mo 2 — l)<l>xx—<f>i/u~ 0**+~2 4>xt — 0 ( 1 ) 

«0 CbQ 

where a a is the velocity of sound in the free stream and t 
denotes time. 

Tke assumptions underlying the derivation of equation 
(1) have been stated in numerous places but are not always 
obviously compatible. It is assumed here that tke ratios 

-^r, are small compared to one, where u, v, w are 

Vo Vo Vo 

induced velocity components and V 0 is the velocity of the 
free stream; moreover, 


TABLE I.— LINEARIZED PARTIAL DIFFERENTIAL 
EQUATIONS OF COMPRESSIBLE FLOW 



Steady state 

Two dimensions 

(A) 

(.vv-i)«„-p„=o, .\r 0 >t, (B) 

qi-.v/^)0„+^ v +* c ,-o, -v u <i, (c) 

Three dimensions i 

,*vr 0 >i, (D) 

Ay+«„-0, M 0 -l, (E) 


Unsteady state 

Two dimensions 

' (0 


The Mack n umb er range for which the equations are valid 
cannot be prescribed a priori since induced velocities are 
functions of wing geometry and angle of attack. We can 
say, however, that for certain configurations at small angles 
of attack the equations and their solutions are consistent 
with the assumptions. In particular, three-dimensional 
lifting surfaces with sufficient- sweepback yield solutions of 
this class at M 0 = 1 . The differential equation shows that 
in this case the boundary conditions need only be specified 
along strips in the transverse direction. The surfaces of the 
Mach cones also are normal to the free-stream direction so 
that any disturbance point makes itself felt at all points not 
upstream of it. Since for these lifting surfaces the disturb- 
ances do not become excessive at M 0 — 1, we have a specific 
kind of lateral strip theory that yields formal solutions 
compatible with the assumptions made. 


7—1 

2 


(2u 


U 2 -tV 2 -tW 


=)«■ 


and, finally, the velocity gradients at a given point of the 
flow field are all of similar magnitude. 

Special cases. — The particular forms of equation (1) to be 
considered are given in table I. In the steady-state equa- 
tions the original independent variables are retained; the 
two-dimensional unsteady-state equation has, however, new 
variables defined by the relations 

x'=x—aoMot, z'=z, i'=a4 (2) 

and consequently the boundary conditions for any particular 
example will be subject to the same transformation. In all 
the equations the constraints imposed by the linearization 
permit, for problems in wing theory, the boundary conditions 
to be specified in the plane of the wing. This plane shall 
arbitrarily be taken to be 2 = 0 . 


BOUNDARY CONDITIONS 

Steady state. —The boundary conditions are given in the 
2=0 plane and in the case of two-dimensional theory the 
wing is assumed to extend infinitely, parallel to the y axis. 
As a convenience in notation, two subscripts will be intro- 
duced: the first, u, denotes conditions on the upper surface 
of the wing, that is, the limit of the function as z approaches 
zero through positive values; the second, l, denotes conditions 
on the lower surface of the wing, that is, the limit of the 
function as z approaches zero through negative values. 

Four types of boundary conditions arise in practice : 

1. Symmetrical nonlifting wing (boundary-value problem 
of the first kind). — The conditions w u =w t = 0 hold over all 
of the xy plane except for the region occupied by the wing. 
On the wing, the relations 2w u = — 2wi=Aw 0 =f(x, y) are 
given, the function/ ( x , y ) being determined by the geometry 
of the configuration. Over all of the xy plane, Av^—u u — 
«i= 0 applies. 
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2. Lifting surface with specified loading (boundary-value 
problem of the first kind). — The condition Au 0 =u u — Ui—0 
holds over the xy plane except for the region occupied by the 
wing. On the wing, the relations 2u„= — %Ui=Au 0 =f ( x , y) 
are given, the function / ( x , y ) being determined by the 
specified loading. Over all of the xy plane, Ato o =0 applies. 

3. Lifting surface with specified camber and angle of attack 
(boundary-value problem of the second kind). — The condi- 
tion Au 0 =0 holds over the xy plane except for the region 
occupied by the wing. On the wing the relation w=f (x, y) 
is given, the function / ( x , y) being determined by the given 
camber, twist, and angle of incidence. Over all of the xy 
plane, Aw o =0 applies. 

4. Symmetrical wing with specified pressure distribution 
(boundary-value problem of the second kind). — The condi- 
tion Aw 0 =0 holds over all of the xy plane except for the region 
occupied by the wing. On the wing the relation u=j ( x , y) 
is given, the function/ (x, y) being determined by the specified 
pressure distribution. Over all of the xy plane, A-u 0 =0 
applies. 

In all cases, induced velocity u is related to pressure 
coefficient C p by the relation ' 



Unsteady state.— The steady-state boundary conditions 
have been given in the most general terms possible. The 
unsteady-state conditions will be limited to a more restricted 
type of problem, namely, cases wherein the airfoil is assumed 
to experience at t = 0 either an abrupt change in angle of 
attack without pitching or starts to travel at the instantane- 
ous velocity V 0 and angle of attack a. In this way the 
transient variation of load distribution and airfoil character- 
istics can be calculated for unit angle-of-attack change. 
Similar methods can treat unit rate of pitching, or deflection 
of aileron, as well as the effects produced when the airfoil 
enters a gust of given structure. The use of solutions of 
such problems in connection with operational methods is well 
known in applied mathematics. Applications of these oper- 
ational methods to aerodynamics have been given by E. T. 
Jones (references 5 and 6) for incompressible fluid theory 
and in reference 7 for supersonic flow. 

If the rectangular coordinate system x', z', t' associated 
with equation (2) is considered to be fixed, the airfoil moves 
in the negative x' direction and the free-stream velocity is 
zero. A simple distortion of the time axis is also introduced 
to simplify the differential equation. Figures 1 (a) and 1 (b) 
aid in the visualization of the problems involved. The air- 
foil section is assumed to lie initially on the x' axis with lead- 
ing edge at the origin and trailing edge at the point k'=c 0 . 
As time progresses the airfoil sweeps across a portion of the 
x't' plane, the leading edge traversing the line x'=M Q t' while 
the trace of the trailing edge is the line x' =c 0 —M 0 t'. The 
region bounded by these lines and the line t' = 0 is that swept 
by the airfoil. The characteristic cones of the differential 
equation cut the x't' plane along lines inclined at ±45° to 
either axis. If the airfoil experiences an angle-of-attack 
change a without pitching, the “area” swept over by the 



(a) Supersonic wing. 

(b) Subsonic wing. 


Fiavus 1.— Diagram for use in determining boundary conditions in two-dimensional 
unsteady motion. 

wing must yield w=—V a a. On the other hand, if the airfoil 
enters a gust of constant vertical velocity w 0 , the region over 
which the modification of w is effective is restricted to the 
region occupied simultaneously by the airfoil and the gust. 
If, for example, the edge of the gust is fixed along the t' axis, 
this axis will form the right-hand boundary of the region 
over which the change in boundary conditions occurs. A 
statement of these boundary conditions may be put in the 
following form: 

1. Lifting surface undergoing abrupt change or starting 
from rest with given velocity. — The condition Auo=0 holds 
over all of the x't' plane except for the region swept across 
by the airfoil. In this latter region, the relation w=f(x',t') 
is given, the function f(x' ,t') being determined by the modi- 
fication in airfoil angle of attack, pitching velocity, aileron 
deflection, or by the gust structure. Over all of the 'a ft' 
plane, Awo=Q applies. 

The expression for pressure coefficient is 



—2 by> 
V 0 M 0 U' 


SOLUTION OF BOUNDARY- VALUE PROBLEMS 

General treatment. — The use of Green's theorem in the 
solution of second-order partial differential equations leads 
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one to the consideration of certain particular solutions of the 
given equations. Because of the physical importance as 
well as the mathematical applicability, attention has been 
centered on the use of a so-called fundamental solution or 
source potential. Thus, in the subsonic case, the potential 
at the point x,z of a unit source located at the point and 
applicable to equation (A), table I, is the logarithmic function 

4>(x, z) = Znr yl = 2 ^g In V(x— £)*+£ 2 (s— f) 2 (3) 


while for equation (G) the potential at x,y,z of the unit 
source at £, 77 + is 


*<**’*) 



(4) 


Here and elsewhere we have (3 2 =[1— iV/ 0 2 [ where the bars 
indicate that absolute values are to be taken. 

The- application of these potential functions to the solu- 
tion of boundary-value problems in subsonic linearized flow 
is well known. Supersonic theory, however, introduces 
added complications when the fundamental solutions are 
considered and, although methods have been established, 
the mathematical techniques are of comparatively recent 
origin. The principal difficulty lies in the integration of 
higher-ordered singularities that appear in the three-dimen- 
sional analysis. Hadamard (reference 4) resolved these 
difficulties and thus avoided the more specialized approach 
of Volterra (reference 8). It would appear, however, that 
a more direct method of derivation stems from Marcel 
Riesz’s use of fractional integrations. (See, in this con- 
nection, references 9 and 10.) The oddness or evenness of 
the number of dimensions still involves considerable differ- 
ences but the final solutions are easily applied. 

In two-dimensional supersonic flow, the potential at the 
point, x,z of a unit source located at the point f,f is defined 
as follows: 


<f>(x,z )= 0 for (x— £) 2 </3 2 (z— f) 2 \ 

1 ? (5) 

<p(ar,z) = — ^ for (x- f) 2 ^ 0 2 (a— 

In three dimensions, the source potential is 

(6) 

at all points for which the radical is real and is zero elsewhere. 

These functions are directly applicable to equations (B) 
and (D) of table I. Equations (E) and (F) are, of course, 
special mathematical cases of equations (A) and (D) for 
which M 0 is 0 and V2, respectively. 

By means of the various source potentials, it is now possible 
to present solutions of the differential equations in terms of 
the prescribed boundary conditions. These conditions are 
assumed to be given in the z=0 plane and subscripts u and l 
shall again denote conditions at 2 = 0 + and s=0— , respec- 
tively. The general solutions appear as follows: 

Equation (A), 

-s#/.' ['» w « (is -I f)- 

(<+-+) (j~j: In J d£ ( 7 ) 

Equation (B), A/ 0 >1 

♦(*.«)— -jj a f » r s>° 

-+iir £***•<» cs> 

Equation (C), Af 0 <+ 

(<pu—4>d dr i ( 0 ) 


Equation (D), ilf 0 > 1 


9{x,y, z) = -±-\jJ [(^) f=0 ** (| aj 


In the last equation, the range of integration is confined to 
that portion of the 2=0 plane that lies within the Mach 
forecone of the point x, y, s, that is, within the half-portion 
of the right circular cone 

( X _ 02 _^[ ( 2 / _^) 2 + ( 2 - r) 2 ] = O 

lying upstream of the point x, y, z. The semivertex angle /j 
of this cone is the Mach angle and is given by the relation 

ji=arc sin -=^=arc cot jS 

The symbol | was introduced by Hadamard and denotes 
the “finite part” of the integral. As in the case of Cauchy’s 
principal value, an improper integral is reduced by a pre- 
scribed technique to a finite and unique value. By definition 
(see also reference 11), 


j A{x)dx __ f*aA(x)—A(x 0 ) , 2+(r 0 ) 

|J. (* 0 -x ) 3 ' 2 (*o -a ) I/2 U ^ 

In the two-dimensional supersonic case, the solution for 
the velocity potential is expressed as the integral of a dis- 
tribution of source potentials. In all other cases, both the 
source potential and its derivative appear in the integrand, 
this latter expression being identified with the doublet 
potential. 

Non lifting case (symmetrical wing, boundary-value prob- 
lem of the first kind). — Equations (7), (8), (9), and (10' are 
applicable directly to the calculation of the potential func- 
tion corresponding to a symmetrical wing. The relation 

dz 

<t> a =4n follows from the condition u u =u t . Moreover, if 
denotes the local slope of the upper surface of the wing, 
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1 50 „ 1 bq>i = /dz u \ 

Vo 5z 5s \dx ) 

and the solutions of the various equations are expressed in 
terms of source distributions alone. 

For example, equation (10) becomes 

V« C C dz„ d£dy 

4>(x,y,e)- ~J r J -fc ^[C»— (12) 

where the finite part sign is dropped since the integral is 
proper. This equation was given by Puckett in reference 3. 
The pressure coefficient on the surface of the wing is 



lifting case (boundary- value problem of the first kind). — 
From the condition w u =w t we have A«j o =0 and the inte- 
grands in equations (7), (9), and (10) are expressed solely in 
te rms of doublet distributions, while equation (8) yields the 
result that conditions on either side of the wing have no 
effect on the other side. 

T aking equation (10) as an example, the solution under 
the prescribed conditions is 


, , 1 I C C , , , \ P 2 zd£dr, 

4>{x,y, «) — 5 F|JJ {{x -t.y-p[( : y-. n y +z *\yn 

(13) 

where 

0a — 0 i= m h« 0 (£> v) d£ 

A more direct evaluation of perturbation velocity u can be 
obtained from the alternate expression 


1 { 

' f Au 0 (£,lj) fizd^d-q 

2tt J 

J {{x-ty-my-nY+z 2 ])*' 2 


Similar expressions exist for equations (7) and (9). 

Lifting case (boundary-value problem of the second 
kind).— This type of boundary condition cannot be solved 
directly by means of the formulas which have been presented 
but resolves always into the required solution of an integral 
equation. In three-dimensional subsonic wing theory, the 
method of solution depends usually on sorqe modification of 
Prandtl’s lifting-line theory although, more recently, lifting- 
surface theories by Falkner (reference 12) and Cohen (refer- 
ence 13) have been applied successfully. 

In the case of three-dimensional supersonic-wing theory, 
sources, s inks , and doublets, have been utilized in two ways 
in the solution of lifting-surface problems. The first of these 
methods was given by Evvard (references 14 and 15) and is 
particularly powerful when one of the leading edges of the 
wing is of the supersonic type, that is, when the velocity 
component of the free stream normal to the edge is greater 
than the speed of sound. A second method of solution was 
presented in reference 16 for the important case of wings 
with subsonic leading edges, provided the flow field about 
the wing is of the conical type introduced by Busemann 


(reference 17). The essential feature of this method is the 
use of differential lifting elements carrying a constant load 
and designed for use in conical flow fields. The solution 
consists of determining the distribution of loading over these 
elements so that the resultant induced vertical velocity at 
any point on the lifting surface satisfies the local boundary 
condition. When approached from this standpoint, the 
problem again requires the solution of an integral equation 
but the equation is of the form 

•»-r*s£ <i5> 

and is thus well known from low-speed airfoil-section theory. 
Inversions of this equation have been provided by Allen 
(reference 18) and von Mises and Friedrichs (reference 19). 



Figure 2.— Lifting-surface element carrying constant load. 


Figure 2 shows the elemental lifting surface to be used. 
The sides of the element extend back from the tip of the 
Mach cone, making angles S and 5+ AS with the positive x 
axis or free-stream direction. The vertical velocity induced 
at the point x, y, 0 by the element will be a function of 5, 
V 

AS, and “ or, changing the notation, 6 , AO, and w where 

8=1 3 tan 5 

6+Ad=p tan (5+A5) 

w=j8 — 
r x 

Denoting the gradient of this vertical velocity with respect 
to 8 by it follows that 

v){8+A8,u)—w(8 ) &>) 

50 Aff-)0 A8 

%■ 

where w(0,ea) and w(0+A0,w) are the velocities induced by 
right-triangular lifting surfaces with constant loading and 
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with vertex angles equal to 5 and S-i-AS respectively. The 
velocities induced by the constantly loaded surface are 
determined directly from equation (13). The results of 
these calculations yield the following expressions: 

For o><C0 


and for w>0 


&t£g=o = Cj 3 r* 

50 2ttVoJ-i 


dwz =0 = Cj 3 r°> 

50 2irT a Ji 


VI— 


c=)i(«i— 6) 


2 da i 




a i 


a^ai—oy 


. dui 


(16a) 


(16b) 


the x axis or free-stream direction, the sides of the triangle 
make the angles S 0 and 5i so that the total vertex angle is 
£>o+S;= 2A. The quantities 0 O > $ lf and a are also introduced 
where 

6fj=P tan S 0 


0i=/3 tan Si 


gj=j3 


y 

X 


Mo 


where 

n-Yl^P 

2 a 

A® • 

The term — is the load coefficient- and is equal to the differ- 
2 

ence between pressure coefficients on the lower and upper 
surface of the wing, 

A P Pi Vu /-i p 

g — g — L P{ 

UTonlifting case (boundary-value problem of the second 
kind). — In the previously discussed lifting case, the induced 
vertical-velocity field of a constantly loaded element was 
calculated. An analogous type of element can also be 
developed for use in the determination of nonlifting wings 
with prescribed pressure distributions. It is apparent that 
differential expressions similar to equations (16a) and (16b) 
must be derived which establish the induced field of the a; 
component of perturbation velocity for a conical-flow ele- 
ment with constant vertical velocity. From equation (12), 
the- following expressions result : 


For cd<0 

bC v 

be 

2X 

f" GJldctft 

J-l (u t — 0) 2 V 1 — (Oi 

(17a) 

and for «>0 

5(7„ 

d0 

_2X 

fa . 

Glides 

(«! — 0) 2 V1 — Cdi 2 

(17b) 


where pressure coefficient C P and surface slope X are 




Figcbe 3.— Yawed triangular flat-plate lifting surface with subsonic leading edges. 


The vertical induced velocity at any point on the wing 
can be found from equations (16a) and (16b), provided the 
distribution of the loading factor C is known. Setting 
(7= (7(0), the downwash is given by the expression: 


Wz 


2irVr 





V l — ar 



(18) 


Since the lifting surface is flat, the function (7(0) must he 
found such that w z=Q is independent of a for — 0 L <ca<Jo- 
The integral equation can be greatly simplified by inte- 
grating the variation by parts and then taking the partial 
derivative of both sides of the equation with respect to a. 
In this manner, equation (18) reduces to 


The application of equations (17a) and (17b) to the 
determination of a thickness distribution supporting a given 
pressure distribution consists of determining X as a function 
of 0 such that the desired pressure results. The essential 
simplification of the method is brought about by the use of 
elements that lead to single integral equations of standard 
form. 

APPLICATIONS 


n _ A yi— m 2 p o C{8)de , i p. C{e)de 

5ca a J—9 ! a — 0 y 1 — a — 0 


which becomes 


<»> 

The solution of equation (20) can be written 


YAWED TRIANGULAR LIFTING SURFACE 

Consider a yawed, triangular flat plate with subsonic 
leading edges such as is indicated in figure 3. Relative to 




A+T, 


i=l (0+Xi) 


0 V(0"T0l) (S—So) 




( 21 ) 
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and, if the integrated loading of the wing is finite, 


0(8) = 


Ae+B 

V(0+0i)(0-0o) 


( 22 ) 


Substitution of equation (22) into equation (18) yields the 
two relations 


and 


where 


and 


Wz=0 " = 2T^o ft) + BH ^> ft)] (23) 

Wz=°=2po [“^(ft, ft) +BH 2 (d u e 0 )] (24) 

H l (d u 8 Q ) = j i —j==^ == — = (25) 

H i (e 1> e 0 )=j i == (ai+flb) -- (26) 


Equations (25) and (26) may be integrated by the standard 
methods for elliptic integrals and, after substituting into 
equations (23) and (24) and solving for A and B, one gets 


and 


where 



(27) 

(28) 


^ . l + ftft-Va-Qflq-ft 11 ) 
ft+ft 


(29) 


and E' is t he com plete elliptic integral of the second kind with 
modulus Vl — G 2 - 

The load distribution over the wing can now be calculated 
from equations (22), (27), and (28). It follows that 


Ap__2C(6)_2a / 2G T(ft-ft)«+2ftftl 

2 tv wMm, Lv(ft+e).(ft-0)j (30) 

Typical load distributions over a yawed wing are shown in 
figure 4 for /3 tan 8i=0.6 and for j3 tan 5 a equal to 0, 0.3, 0.6, 
and 0.9. 

Integration of equation (30) over the surface of the 
triangular wing determines for lift coefficient the expression 


n 2 ax /fftanA 

G l=W cos A Y p (31) 

where A is sideslip angle and 2A is the angle between the 
leading edges. Equation (31) was derived for wings with 
subsonic leading edges and supersonic trail ing edge and 
consequently is valid only for cases for which 


ju+A < 90 
A+AO 
A — A > 0 


Mo 



- 1.0 -8 -.6 -. 4 . -.2 O .2 .4 .6 .8 tO 

3 ton 3 


Floras 4.— Angle-of-attack load distribution over yawed triangular plan form, p tan 5] —0.0. 


DOWNWASH BEHIND TRIANGULAR WING IN SUPERSONIC FLOW 

The second application will show how doublet distributions 
may be employed in the calculation of downwash in the wake 
of an unyawed triangular wing with subsonic leading edges. 
The expression for the velocity potential will be given in 
general form, but in order to avoid detailed analysis the value 
of downwash is determined only along the center line in the 
plane of the wing. 

\m 0 



A plan view of the wing and wake is shown in figure 5. 
The load distribution over the wing is found from equation 
(30), after setting 0 o =ft , to be 


( 32 ) 
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4 P- 46 0 2 ax. 

2 E'fjafr-t 3Y • ; 

where E' 0 has the modulus ^1—8*. Moreover, 

<p— j u dx 

so that the jump in potential 0,— 0 r on the surface of the 
wing is 

V^ 2 -) SV (34) 

and for points in the wake is 


Region. B : 
where 


w 2E X t 2 pi K—E 

Wo ‘S’E'o ^irE'o Jo h-\-8o 


(38) 


w 0 induced vertical velocity on the wing 

K complete elliptic integral of the first kind 

K X ,K« complete elliptic integral of the first kind with 
moduli A*i, k lt respectively. 

E complete elliptic integral of the second kind 

E X ,E 2 complete elliptic integral of the second kind with 
moduli 1 * 2 , respectively 


0a“ 


2 aVp 

o/3 


-4>i=gr-; VW-i8V 2 


(35) 


l*i 


Cq 6 q 
X — Co 


where c 0 is the root chord of the wing. Since the wing and 
its wake form a discontinuity surface for the perturbation 
velocity potential and since for all points on this surface 

50, 502 

5s 5s 

it follows from equation (10) that the velocity potential at an 
arbitrary point x,y,z is given by the relation 


, X—Cp 

2 c 0 8 0 

V) 

Figure 6 shows the variation of — along the x axis for var- 

ious values of the parameter 9 0 =l 3 tan 5„. The asymptotic 
values at x= c° are also indicated and can be shown to agree 
with the values of downwash at infinity for a wing with the 
same span load distribution in incompressible flow. The 
discontinuity in downwash at the trailing edge is a character- 
istic property of supersonic-type trailing edges. 


0(z, y, s) j/J (0.-0J (|: £) r _ o dtdr, (36) 

where r is that portion of the wing and wake forward of the 
Mach forecone from the point x, y, s and 0,-02 is given by 
equations (34) and (35). 

The value of the downwash aft of the wing and along the 
x axis will he calculated from equation (36), thus 

V~Q 


is to be determined. In carrying out these calculations, it is 
necessary to consider two segments of the x axis behind the 
wing: 

Region A extends from the trailing edge to the point where 
the trailing Mach cones from the tips of the wing intersect 
the x axis and thus includes values of x satisfying the in- 
equality 

Co^a;^eo(l+0o) 


Region B includes values of x for which 



Figcbe 6 —Variation of downwash on x axis behind a triangular wing plotted as a function 
of distance in chord lengths. 


Denoting downwash at this point by w,, Lagerstrom and 
Graham (reference 20) have shown that 


w t E' 0 —d 0 
w 0 E'o 


(39) 


Co(l-Mo) ^x 

The final expressions for downwash in the two regions are 
found, after some manipulation, to be 

Region A: w n E 2 — (1— h 2 )K> , E' 0 — 8 0 
■w 0 ~ irE'Jc, -r E'o 

2 K—E „ 

rtf'a Jo k\l+8jc) dk (3/) 


A more detailed development of the results shown in figure 
6 has been given in reference 21. 

TWO-DIMENSIONAL UNSTEADY LIFT PROBLEMS 

It has already been pointed out that, in the case of unsteady 
motion in a two-dimensional compressible-flow field, the 
linearized partial-differential equation for the perturbation 
velocity potential can be transformed into the same form 
that has been considered in solving steady-state problems 
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in supersonic wing theory. This immediately indicates the 
possibility that for certain types of boundary-value problems 
in the unsteady case an analogy can be established with three- 
dimensional lif ting-surf ace problems. 

As an example, consider an airfoil that has been flying 
at supersonic speed and then, experiences at t'—O an abrupt 
angle-of-attack change without pitching. Since the angle- 
of-attack change is assumed to take place at f'=0, it can be 
assumed that previous to this time the induced velocities of 
the wing are zero and only subsequent perturbations are 
to be calculated. Throughout the swept area in the x't' plane 
(fig. 1) the vertical induced velocity w is constant and equal 
to —TV*- Elsewhere in the s'=0 plane there is no dis- 
continuity in the value of pressure, that is, is continuous 
at z'= 0. 

Suppose now that the area is a wing plan form and that 
the free stream is directed along the t' axis. The character- 
istic cones of the unsteady problem become the Mach cones 
of the steady-state problem, and the Mach number of the 
free stream is V2 since the characteristic lines in the figure 
are inclined 45° to the axes. Moreover, the induced ver- 
tical velocity is ^ and the perturbation velocity in the free- 

stream direction is A correspondence can thus be estab- 
lished between the unsteady problem and a three-dimensional 
lifting-surface problem. 

As outlined, the boundary-value problem is one of the 
second kind, that is, w is specified on the wing and 
Au=Aw=0 off the wing. In this particular example, how- 
ever, the edges of the wing are of the supersonic type and 
no interaction exists between the two surfaces of the lifting 
plate so that pressures on either side can be. calculated by 
the methods used for symmetrical nonlifting wings. Thus, 
from equation (12), for s'>0. 

*') -gy jvy-i’ (4o > 

and for all z' 


x', z') = <p{t' , x'-z ') 


Ap 


The expressions for the indicial load coefficient are as 
follows: 


Region A (between lines x'=t', f'= 0, and x'=Co—M 0 t') 

Ap 4a 

g 


(41a) 


Region B (between lines x' = —t', x'=t', and x'=e 0 — M 0 t') 


Ap_ 4a 


Mdx'+t' | VMq 2 — 1 


2 ViW-l [it arc cos x' +M 0 t' ^ tM 0 

^+arcsin|r^J (41b) 

Region C (between lines — and x' =—t' 

(41c) 


The growth of -£■ with time, as obtained from equations 

(41a), (41b), and (41c), is shown in the portion of figure 7 
designated “supersonic”. At £'= 0 the loading jumps to the 

value jj- and is constant along the entire chord. This value 

persists throughout the previously denoted Region A and 
thus, with advancing time, moves rearward along the chord, 

leaving the trailing edge at Over the forward 

portion of the chord the familiar Aclieret type of steady- 
state loading becomes effective, spreading back from the 
leading edge and occupying the entire chord length after 


t' 


An 


Previous to t f 


Co 


a transition region be- 


M l JL 1 U Y A \J UO UW C 71 ,r | -t 

0—1 Mo+1 

tween the two types of constant loading exists, and subse- 
quent to this time this transition region moves aft and 

c 0 


leaves the trailing edge at 


qcC 


Ve. 

qcC 




Figure 7.— Pressure distributions on wing undergoing sudden onglo-of-attack change at i'*> 0. 

For purposes of comparison, the growth with time of the 
angle-of-attack indicial load coefficient for subsonic flight is 
also shown in the part of figure 7 entitled “subsonic.” Since 
in this case the lifting-surface analogue involves subsonic 
leading and trading edges, the analysis requires the solution 
of a boundary-value problem of the second kind. The 
method of Evvard (reference 15) was used to obtain the 
results shown. It is to be noted that the expression 

Ap 4a 

~q~Mo 

holds at t'= 0 for all values of Mach number. 

Figure 8 shows the variation of the indicial lift function 
C L (0 defined by the relation 


a 


, ( o =- 1 r 

' C 0 aJ 0 


Ap 


dx 


as a function of Mach number and half-chords s 


( 42 ) 
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t'+Afao-t'Mo* 


SMgt' 



Figure 8. — Indicia! lift-curve slope for Mach numbers between 0 and 1.4 shown to time 
required to travel 12 half-chord lengths. 

traveled. The curve at M 0 =0 was first studied by 
Wagner (reference 22) and R. T. Jones (reference 6). Since 
4 

the starting value is C La (t') must initially be infinite. 

Immediately afterward, however, it assumes the value ir and 
then rises to the asymptotic value of 2-ir. At a Mach number 
•of 0.4 the starting value of C La 0 1 ') is 10 followed by a de- 
decrease for the time required to travel approximately one- 
half chord length and finally a steady rise takes place to the 

asymptotic value ~-.f i’ At M 0 =0.8 the behavior is 

similar. The dashed portions of the curves were dete r mi n ed 
from the known variations of the functions and were not 
calculated explicitly. The asymptotic values of C r ia con- 
sistent with the Prandtl-Glauert correction become so high, 
however, with increasing Mach number that the assumptions 
of small perturbation theory are undoubtedly invalidated 
near M 0 =l for sufficiently large values of s. The initial 
portions of the subsonic curves shown in the figure are, 
however, valid results of the theory. The nature of the 
indi c ia! lift function is somewhat different at supersonic 
Mach numbers in that the beginning portions of the curves 
are flat. The curves rise afterwards, however, in a finite 
time to their steady-state value. From equations (41a), 
(41b), and (41c) the expressions for C La (f') are easily cal- 
culated for Mo Si and are as follows: 


arc cos 


First time interval 0<C< 

4 




Mo 


Second time interval ^-rxr<f / < i C ° 


(43a) 


1+Mo ^Mo-1 

. c 0 —M 0 t'' 
sm - — 


[sfc (l +aic 811 c -^~r £ ) 


VW=T 


Cq 


(43b) 


Third time interval 


Co 

M 0 -l 


<t' 




Vm 0 2 -i 


(43c) 


Some of the above results, along with further developments 
involving the entrance at supersonic speed of an unrestrained 
airfoil into a gust, have been given in reference 7. 

SWEPT-BACK WINGS AT ATo=l 

Consider now the special form of the basic differential 
equation for the case M 0 =l. As given in table I, equation 
(E), the velocity potential satisfies Laplace’s equation 

in two dim ensions. The boundary conditions need, there- 
fore, to be given along strips normal to the free-stream direc- 
tion. Equation (7) expresses the solution of the equation in 
terms of two-dimensional sources, sinks, and doublets where, 
in these variables, 

In rj_=ln-yl(y—o) 2 +(z—fi 2 

The proposed problem is the determi n ation of the angle- 
of-attack load distribution over a swept-back lifting plate, 
the leading edges will be assumed straight lines while the 
trailing edge will, for the time being, be left arbitrary. The 
nature of the wing is thus indicated somewhat arbitrarily in 
figure 9 (a). 

Denoting the semivertex angle by S 0 so that the equa- 
tions of the leading edges are 

y=±r tan 5o=±mj: (44) 

it follows from equation (7) that since 

S<P„ b (jil 

dz ~~ bz 

the velocity potential is given by the relation 


where 


, 2 C™ A 4> 0 (x,v) J._ 

*{x,y,z) = 2?r j-mc (y-ny+¥ dv 


Atp 0 (x,y)= \ Ai? 0 (j;,i7,0)di7 


(45) 


It is then possible, after integrating equation (45) by parts 
and imposing the condition that A <f> a (x,y)=0 at -q—±mx, to 
calculate the derivative of 4> with respect to z and thus obtain 
for w a , induced vertical velocity on the wing, the expression 


— 1 f mx Ai’ofo n) J 
Wo=-nY - m 
Zn j—TBx y v 


(46) 
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(a) Plan form with arbitrary trailing edge. 

(b) Plan form satisfying Kutta condition. 
Figtjse 9.— Swept-back wings for analysis at Afo=l. 


This integral equation is to be solved for Afi 0 , the velocity 
w 0 being assumed constant on the wing and equal to — V 0 a 
where a is angle of attack. The load distribution can then 
be calculated from this solution by means of the relation 


t\p _2 l 
i 


A 


d<fo_ 2 d 
lx ~F 0 hr 



(47) 


The remainder of the analysis can best be divided into two 
parts: The first case treating values of x between 0 and Co, 
the second dealing with the remaining values of x on the 
wing. 

Case I: — Since the leading edges of the wing 

are of subsonic type, singularities in pressure occur at these 
edges so that the required solution of equation (46) is of the 
form 


Av 0 (x,y) = 


A+By 

■yjrtftf—y 2 


(48)' 


Substitution into equation (46) and use of the fact that 
A<j) 0 (x, y ) is an even function of y leads to the result A= 0, 
B= — 2w 0 . Hence 

(49) 


and 


Ap _ —4 w 0 m 2 x _ 4am a r 
2 FoVmV-y 2 ^!mV—y 2 


(60) 


Case II: c 0 =£- — Let the equation of the trailing edge be 


y—a(x) or x=a*(y) 

Using in equation (46) the fact that 

A»o (*, y) = —Av a (x,—y) 
the expression for w 0 becomes 

— 1 C a(x) 2yAv 0 d-rj 1 C m 2rjAv 0 dr) 
^=27 jo 


(51) 


(52) 


If, on the surface of the wing, Afo{z,t\) is known, then, in 
the wake, the discontinuity in the velocity potential is 
A 0 o[&*Oj),»?] since no contribution to the jump in potential is 
made past the trailing edge. It follows that if on the wing 


Av 0 (x,ti) =f(x,y 2 ) (53) 

then, in the wake 

(53) 

M V, So* + y [tl , w _ (54) 

Substituting equations (53) and (54) into equation (52) and 
introducing the Kutta condition that loading at the trailing 
edge is zero, one gets the modified integral equation 




(55) 


where the variables a, <?i now replace y 2 , respectively. 
The function 


f(x,<r) = 2w Q 


V 


or— a 2 (x) 
m 2 x 2 —<r 


(56) 


satisfies equation (55) and it remains to determine a(x) so 
that pressure is zero at the trailing edge. But from equation 
(53) 

iy, 

and thus 

(57) 

where k'=-^~ and E{k,ip), are incomplete elliptic 

integrals with modulus k=^l—k n and argument ^=arc 

• i r, ¥~ 

sm 7 - -»/ l—-~n 

k V mx 


AM%, y) =— 2w 0 VmV— v 2 
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At the trailing edge y=a(x) and 



where the elliptic integrals are now complete. If load co- 
efficient is set equal to zero, the differential equation 


da m?x E 
dx~ a K 


(58) 


Ox — SVo A4>(T.E.,y)dy (59) 

where A 4> ( T . E., y) is the jump in potential at the trailing 
edge and thus equal to the circulation function r (y). The 
following results are obtained: 


r (y) — 2mV a c 0 a, 0 £y£A 

T(y)=2V a ab [E(k , *) - (1 -P) FQc^A^y^b (60) 


/ d“ 

where &=-y follows. 

The integration of equation (58) leads to the shape of the 
trailing edge for which the Kutta condition is satisfied. 
Figure 9 (b) shows the plan form of the wing. It can be 
shown that the slope of the extended trailing edge approaches 
the slope of the leading edge. 



Figure 10. — Load distribution, over swept-back plan form at A/o=I. 


In figure 10 the load distribution is shown at three span- 
wise stations for the case when the wing is cut off along a 
line normal to the free-stream direction. Over the center 
section of the wing the Ackeret type of distribution exists. 
The remaining sections have discontinuities in slope of the 
loading at the point where the chord line is cut by the Mach 
cone arising at the trailing edge of the root chord. This 
behavior of the loading has been noted elsewhere for swept- 
back wings at higher Mach numbers. (See, e. g., reference 
23.) 

Lift coefficient C L of the wing is given by the expression 


c "-hfjyS^t dx 


where S' and b are, respectively, area and semispan of the 
wing and the first integral extends from the leading edge to 
the trailing edge. This equation may be rewritten as 


where k 




and A is the lateral distance to the in- 


board tip of the wing. (See fig. 9 (b).) In figure 11 the value 
of 



0 -S .4 .6 .8 1.0 I.S 1.4 

y/co fan 6 0 


Figure ll — Spanwise distribution of circulation for swept-back wing at Afo=l, 6=1.325 co 
tan So and .112=4.57 tan Si. 


r y 

Yf r r is plotted as a function of — ~ — 5 - for a wing with 

VaCaa tan 5 0 r c 0 tan S 0 b 

semispan 6=1.325 c 0 tan S 0 and aspect ratio AI2=4.57 tan 

8 0 . Results of the integration of equation (59) are shown in 


figure 12 where 


tan5 0 


is plotted as a function of 


AB 
tan 5 0 



AR/fan 60 

Figure 12.— Lift-curve slope as a function of aspect ratio for swept-back wing at .V/a=l. 

The methods presented here can be applied to the case of 
the swept-back wing with tips cut off parallel to the free 
stream. In this case a Mach cone originates not only at 
the trailing edge of the root chord but also at the intersection 
of the leading and the lateral edges. On the portion of the 
wing downstream of this Mach cone, the load distribution is 
modified so that an abrupt discontinuity exists at the Mach 
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cone and zero loading is effective over this part of- the wing. 
A similar effect on this type of swept-back wing has also 
been noted at higher Mach numbers in reference 23. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Dec,. 22, 1948. 
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